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o , 

, We show that, with reasonable hypotheses leading to a simple modeling, a link can be obtained 

■ from experiments on the axial low frequency collective modes between the molecular scattering 

length um and the energy parameter ^ = 1 + /3 of the gas at the unitarity limit. We also point out 
that, in order to reach the range where the features of the Bose limit can be clearly seen, experiments 
have to go to more dilute situations than have been achieved presently. 

in , 

Collective excitations of an ultracold dense ^Li Fermi gas, most likely in the superfluid regime, have been studied 
experimentally for the first time very recently [1,2]. In particular Bartenstein et al [2] have covered in part the 
crossover range going from the unitarity limit to the BEC regime domain [2]. This is the region to which we will 
O restrict the scope of this paper. In this domain the scattering length a is positive. Physically this is the range where, 
^ for a dilute gas, molecules form out of the fermions. At unitarity where a = oo these molecules are infinitely large 
' and get increasingly small when a decreases. In this last limit these molecules are expected to behave as bosons, and 
indeed the Bose-Einstein condensation of these molecules has by now been observed by several groups [3]. On the 
^ other hand this will obviously no longer be true at unitarity. Hence studying the crossover where the internal stucture 
'■^ , of the molecules will manifest itself progressively is a very interesting question. 

' In this paper we consider the detailed information which can be extracted from the experimental results on the low 
O \ frequency collective modes in this crossover [4]. We show that, in order to see the specific features of the Bose limit, 
experiments have to go to more dilute situations than what has been done up to now. With reasonable hypotheses we 
find that a link can be obtained from experiments between the molecular scattering length aj\/ and the energy of the 
^-H , gas at the unitarity. Hence experiments could provide qm together with the first corrective term in the energy to this 
^ ■ dilute limit, which is of very high theoretical interest. Moreover we show that, even in the range of dilution which has 
^2 , already essentially been reached, precise experiments can already extract these informations from a detailed analysis 
' of the data. More specifically we will consider in this paper only the quasi one-dimensional situations found for very 
J — ] elongated trap clouds [1,2]. However our approach could be extended to other less anisotropic situations. We restrict 
• also ourselves to the case of atoms being only in two hyperfine states, with equal number of atoms in each hyperfine 
] state. For practical matters we will consider only ®Li, on which all mode experiments have been performed to date, 
i In our range of interest we expect superfluid hydrodynamics to be valid because the energy scales set by the 
molecular binding energy or the pair breaking energy can easily be large, compared to the mode frequencies. In 
a related way experiments have already reached temperatures low compared to these scales, so that there are few 
y ' elementary excitations, dissipation can be neglected and we can use perfect fluid hydrodynamics. Indeed dissipation 
' , is found to be quite small [2] experimentally. Hence this domain is expected to be much simpler than the negative 
(— ( a region, where the pair breaking energy will become smaller as |a| decreases. In particular this will make damping 
Q higher, as it is indeed found experimentally. Therefore the validity of reactive superfluid hydrodynamics is much more 
O ' uncertain when one goes in the negative a domain. 

' The mode frequency is determined from the dependence [5] of the chemical potential ^i{n) on the total atomic 
density n. More precisely it is related to the derivative of ^(n) with respect to n so it conveys detailed and sensitive 
information on /i(n) and on the energy of the gas. It is conversely possible [6] to invert experimental data on the 
mode frequency to obtain the equation of state of the gas. This requires only that /.t(n) is known in some limiting 
case, from which one then go away by an iterative procedure, making use of the experimental knowledge of the mode 
frequency as a function of particle number. The basic ingredient of this method has been shown recently [7] to have 
an accuracy of order lO^''. However presently the number and the accuracy of the data is not enough to carry out the 
above program with a sensible precision. Nevertheless we will show that, with reasonable hypotheses, it is possible to 
analyze the data and extract quite important physical quantities. 

We will assume that the (homogeneous) gas can be completely characterized when its density n and the scattering 
length a are known. This last quantity is controlled experimentally through the applied magnetic field. This hypothesis 
is far from obvious since one might expect that, at least in the vicinity of the Feshbach resonance, more parameters are 
needed to fully characterize the scattering properties. However, at least in the case where the width of the resonance 
is quite large compared [8] to the Fermi energy, this seems to be valid. This situation is found in ^Li and in *°K. Then 
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for homogeneity reasons the T = chemical potential of an homogeneous gas with density n takes the general form: 



H{n) = Epfil/kpa) (1) 

where Ep = fi^kp/2m is the Fermi energy and kp = Sw^n. This implies in particular at resonance = a universal 
behaviour since the only energy scale is the Fermi energy. In this case fj,{n) = £,Ef where ^ = 1 + /3 = /(O). The 
energy per atom of the gas itself is given by {3/5)£,Ef- 

Actually in an harmonic trap the experimental determination of the mode frequency as a function of scattering 
length allows to obtain only within a multiplicative constant. This is for example well known for the Bose case 
where only the linear dependence [9] of /i(n) = g'n on n is necessary to obtain the theoretical result for the frequency, 
but the constant g = 4:TTh^a/m drops out, so the mode frequency docs not allow to reach the scattering length a. 
Similarly in the unitary case [10] the frequency docs not depend on ^ and one has rather to find ^ from measurements 
of the interaction energy of the gas [11], as it is obtained from expansion experiments. However since there is only a 
single multiplicative constant, one could hope to find from the mode frequencies the scattering length a as soon as ^ 
is known, or the reverse. Naturally in order to do so, one needs to make a connection between the regions with small 
a (in the Bose regime) and with large a (near the unitarity limit). We will make this connection by taking a quite 
reasonable modeling between these two limiting regions. 

In the vicinity of the unitarity limit w 0, it is natural [12,13] to proceed to an expansion of f{y) since y — l/kpa 
is small. Hence we write: 

m^^^-Sy (2) 

Near the other limiting case, namely the Bose limit, it is also natural to proceed to an expansion [10,14] in powers of 
the density. The standard expression [10] of the chemical potential in the Bose limit is /zm = 2/i(n) = gunu in terms 
of the molecular density Um and the coupling constant is given by qm = ^irTi^aM /rriM in terms of the molecular 
mass JTiM = 2m and the molecule-molecule scattering length ajv/- The next order correction in the standard Bose gas 
expansion is [14] the Lee, Huang and Yang (LHY) term {^2/^y/TT)gMnM{nMa\iY^^ ■ Hence we note that the expansion 
which comes out naturally in this theoretical framework is in powers of n^^^. On the other hand, as we have seen 
above, it is rather the variable kpa ^ v}!'^ which comes out in studying the BEC-BCS crossover, corresponding 
physically to the fact that the internal fermionic structure of the Bose molecule under consideration is now coming 
in. It is unclear whether n^/^ terms will be present. At least it is an open theoretical question to find the specific 
form of the low density expansion, and in particular whether only n^/^ or v}!'^ or both are coming in. In particular 
experiments on collective modes could contribute to clear up this question. In the present case, since our purpose is 
to make the connection to the Fermi gas regime, we want naturally to keep the n^l^ variable. On the other hand we 
will not include for simplicity v}!"^ corrective terms, whose presence is unclear. They are anyway of higher order than 
v}^^ . In this way we are led to write the low density expansion of our function /(y) as 

f{y)^{Aly)[\ + cly) (3) 

This corresponds to the expansion /x(n) ~ (1/2)5m^m(1 + c(37r^)^/"^n^/^a), provided we set A = {l/iTT)aM / o- 
Naturally parameter c is physically important since it is responsible of the departure of the mode frequency from its 
ideal Bose gas value in the vicinity of the Bose limit. Similarly parameter S is responsible [13] for the linear behaviour 
[2] of the mode frequency in the vicinity of unitarity. 

We will now connect smoothly the unitarity domain and the Bose domain by making use of a Fade approximation 
[15] for /(y) with y running in the domain [0,oo[. This kind of approximation is known to be often quite good in 
reproducing a function because it has the potentiality of imitating the analytical structure of the function by the poles 
it may introduce. Specifically it amounts to approximate the function of interest by a rational fraction. In the present 
case it happens that we can make in a natural way such an approximation without introducing more parameters than 
the four ones we have already discussed. Indeed we will take the following Fade approximation: 

f{y) = A y^^l (4) 

y^ + qiy + qo 

Requiring this form to give the proper expansion for small and large y leads to the following expressions for po, qo and 

qv 

^ A + ^c Apo 
Po = ^ ^2_^S 10 = — qi=Po-c (5) 

Naturally behind this approximation we make is hidden the assumption that /x(n) is a smooth function of n when one 
goes from the unitarity to the Bose regime. This means that there is no break in the physical properties during this 
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crossover. Although this is a most common assumption, it is by no means obvious and it is important that it can be 
justified by detailed experiments on the mode frequencies as a function of scattering length. Indeed it could be then 
seen if Eq. 4 gives or not a proper description of these frequencies throughout the whole range of scattering lengths. 

We could naturally treat all of our four parameters as free parameters. However this would make our study 
somewhat heavy. Moreover, in the case of ^Li, there are already rather reliable informations on the parameters ^ 
and S. Indeed experimental [11] and theoretical [16,17] evaluations of ^ converge around a value ^ « 0.45. Whatever 
the final outcome it is not expected to be so far from this value. So we will take ^ = 0.45 as granted. Similarly the 
analysis [13] of the axial mode frequency data [2] leads to 5 w 0.5, so we will take S = 0.5. This leaves us with only 
two parameters A and c. 

Parameter A is essentially the molecular scattering length and it is our central interest. It has been shown theo- 
retically [18] that Om = 0.6a in a nearly resonant situation, which corresponds to ^ = 0.064. It would naturally be 
very interesting to check experimentally this result and we will sec that the mode frequency offers a promising way to 
do it. On the other hand parameter c is the equivalent of the LHY correction in our case and is also clearly of high 
theoretical interest since it should contain contributions due to the fermionic nature of the molecular components. 
On theoretical grounds it is unclear to even decide the sign of c. However if we take into account the analytical result, 
to be given below, for the mode frequency near the Bose limit, its extrapolation toward the unitarity region seems to 
go clearly below the experimental results [2] if c < 0. We will therefore restrict ourselves to the domain c > 0. 




FIG. 1. Reduced axial mode frequency squared as a function of arctan(l/A;Fa) for the model Eq. 4 for c = 0, 2 and 4, and 
values of ttM /ffl as indicated near the various curves. For comparison we have also plotted with the c = case, as the thin line, 
the result of the model of Ref.13 /(«/) = 1/2 — (I/tt) arctan(7ry/2). 



We will now, for consistency, restrict further our parameter domain. First we want f{y) to smoothly interpolate 
between the unitarity limit where /(O) = ^ and the Bose regime where f{y) « A/y, so we do not want any zero of 
/(y) for y > 0. This implies po > 0, i.e. A < 0.4 or UM/a < 3.8 (note that in this case we have also qo > and 
qi > 0). This relation is well satisfied by all the current possibilities [19] for aujo- considered in the literature. Next 
we certainly do not want any pole for y > 0. But even poles for y < correspond to a nearby singularity which 
contradicts somewhat our smoothness requirement. Hence we will require that /(y) has no pole for real y (since we 
have q\ > 0, the real part of the poles location is negative, which puts them away from our y > domain). This means 
qf < 4go- This condition implies that A < 0.3 (i.e. om/o < 2.8) and, in this case, we have the further restriction 
c < (0.73/^)[l + (1 - 2A7Af/^] - 2.7. 

Before discussing our results, let us first briefly indicate how we proceed to get them and give analytical results in 
limiting cases. Since we are interested in the axial mode where the gas behaves as a ID system, we have first [7] to 
calculate, from the relation between the 3D density n(/i) and the chemical potential /i, the relation ni(/i) = dii'n{fi') 
between the effective ID density ni and the chemical potential /z on the long axis of the trap. In our case this calculation 
can be performed analytically from Eq.(2), but since the result is somewhat lengthy we will not give it explicitely 
in the general case. On the other hand this relation becomes much simpler near unitarity where f(y) ~ ^ — Sy. 
Introducing a parameter t it can be written as ni = — {5/8){S/^)t^, with n = — {S/£,)t. In these expressions for 
ni and fj,, we have actually dropped useless multiplicative factors since they do not enter the final result for the mode 
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frequency. Similarly near the Bose limit where /(y) = (^/y)(l + c/y), we obtain ni = + [8/7) ct^ , with fi = + ct'^. 
In the general case analogous parametric relations are obtained. The above parameter is actually given by t = n^/V 
where n is the (3D) gas density on the long axis of the trap, divided by its maximum value found at the center of 
the trap, and k = hm with yM being the maximal value of y = l/kpa found at the center of the trap. Hence the 
maximum value of t is Im — kprnaxO,. 




FIG. 2. Reduced axial mode frequency squared as a function of arctan(l/A;Fa) for the model Eq. 4 with the value of 
clm/ci = 0.6 of Petrov et al and for c = 0, 2, 4 and 8. The filled squares, with error bars, are the data of Bartenstein et al. 



We then calculate the mode frequency by making use of the " a -p modeling" together with a first order correction 
[6,7]. We take the above relation between fii = ni/nimax and p, = fi/fimax (with nimax and fimax being the 
center trap values of ni and fi) and make a best approximation of it by ni = [I — (1 — p,)"/^]^. We then make 
a first order correction to take into account the difference between this model and the actual relation ni(/2). This 
procedure leads [7] to a precision of order 10"'^. At unitarity we have the simple result fii = p-^^"^, equivalent 
to a = 2 and p — 5/2, leading [10] to oj'^/ojI = 12/5. Then our first order correction leads, in the vicinity of 
unitarity, to the analytical [12,13] result ij^juj^ — 12/5 + (256/8757r)(S'/^)(l/fcFmaa:a)- Similarly in the Bose limit 
we have n\ = /I^, equivalent to a = 2 and p — 2. This leads to oj'^ /oj^ — 5/2. Then a perturbative calculation 
gives near this Bose limit uj'^/ujI - 5/2 (3^/^/32)(r(13/3)/r(29/6))c(fci.™aa;a) ~ 0.082c(fcFmai.a). This is to 
be compared with the coefficient one would obtain [14] from the LHY result, which can be rewritten as f{y) = 
{A/y){\ + 64(aM/67ra)3/2/y3/2) ~ {A/y){\ + 0.78(aAf/a)3/2/y3/2). The corresponding result [14] can be written 
uj'^/ujI - 5/2 = (105/256V^)(aAf/a)3/2(fcj.„^^a)3/2 ~ Q.OMiaM / af/'^{kFmaxaf/'^ . Hence the numerical coefficient 
we find in our case, for c ~ 1, is quite similar to the one which comes out in the LHY case. 

We give in Fig.l the results of our numerical calculations for c = 0, 2 and 4, and values of aujo, compatible with 
the restriction given above. In Fig. 2 we take the result cm/o = 0.6 of Petrov et al and give our results for c = 0, 2, 4 
and 8 together with the experimental results. In agreement with our specific choice of S/^ all our curves start with 
the same slope near unitarity. Near the Bose limit our results rise (for c 0) above the Bose limit uP' juJ^ — 2.5. 
The slope for this rise is naturally proportional to c, in agreement with our analytical result. However we notice 
that this rise is always rather small, in contrast to what one might guess at first. Furthermore we see that values 
above the Bose limit are essentially only obtained in the density domain beyond the one which has been reached 
experimentally, so that a positive value for c is not in contradiction with present experimental data. From Fig. 2 we 
see that ajv//a — 0.6 and c > 2. seems compatible with present experimental data. In this respect we notice that when 
we translate a dilution parameter [2] umo^i = 10^"^, which sounds to correspond very much to a dilute situation, 
we find (assuming aM /a = 0.6) l/kpa ~ 1.54 of order unity, and arctan(l/fci?a) ~ 1. is not in the range where the 
typical Bose behaviour occurs. 

Next we see that, starting from the unitarity value w^/w^ = 2.4, the rise toward the Bose value is faster for small c 
and aAf/a. This is easy to understand from Eq. 5 since these small values imply also small values for our parameters 
Po, qo and qi. Accordingly the Bose regime, which from Eq. 4 is found for y ^ po, q^, qi^ is rapidly reached. For larger 
values of c and aM/a, parameters po,qo and qi get correspondingly larger and the Bose regime is reached only for 
larger values of l/kpa. Nevertheless increasing c or um /o. does not lead to the same change in the mode frequency, as 
it can be seen from Fig. 1. Hence by making precise measurements up to l/fci?a — 2, it should be possible to obtain 
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both c and Um /a, and even to check if our approximation Eq. 1 of the equation of state gives an accurate description 

of the data. 

In conclusion, by making use of a simple modeling, we have shown how it is possible, from experiments on the 
axial low frequency collective modes, to find a link between the molecular scattering length om and the energy of the 
gas at unitarity. We have also pointed out the possibility that the specific features of the Bose limit will be seen in 
experiments only by going to more dilute situations than what has been done up to now. 

We are most grateful C. Salomon for stimulating this work and we acknowledge useful discussions with F. Chevy, 
Y. Castin, C. Cohen- Tannoudji and J. Dalibard. 
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